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Abstract 

We consider J/\I/ meson production in ep deep inelastic scattering in the colour singlet 
model using the A>r-factorization QCD approach. We investigate the z-, Q 2 -, p 2 --, y*- and 
H^-dependences of inelastic J/ty production on different forms of the unintegrated gluon dis- 
tribution. The p 2 -- and Q 2 -dependences of the J/ty spin aligment parameter a are presented 
also. We compare the theoretical results with recent experimental data taken by the HI and 
ZEUS collaborations at HERA. It is shown that experimental study of the polarization J/ty 
mesons at Q 2 < 1 GeV 2 is an additional test of BFKL gluon dynamics. 
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1 Introduction 



It is known that from heavy quark and quarkonium production processes one can obtain 
unique information on gluon structure function of the proton because of the dominance of 
the photon-gluon or gluon-gluon fusion subprocess in the framework of QCD [1-4]. Studying 
of gluon distributions at modern collider energy (such as HERA, Tevatron) is important for 
prediction of heavy quark and quarkonium production cross sections at future colliders (LHC, 
THERA). At the energies of HERA and LEP/LHC colliders heavy quark and quarkonium 
production processes are so called semihard processes [5]. In such processes by definition 
the hard scattering scale \i ~ uiq is large compare to the Aqcd parameter but on the other 
hand ji is much less than the total center-of-mass energy: A QCD <C \i <C \fs. The last 
condition implies that the processes occur in small x region: x ~ mg/\/s <C 1, and that the 
cross sections of heavy quark and quarkonium production processes are determined by the 
behavior of gluon distributions in the small x region. 

It is also known that in the small x region the standard parton model (SPM) assump- 
tions about factorization of gluon distribution functions and subprocess cross sections are 
broken because the subprocess cross sections and gluon structure functions depend on a 
gluon transverse momentum [6-8] . So calculations of heavy quark production cross sec- 
tions at HERA, Tevatron, LHC and other collider conditions are necessary to carry out in 
the so called /c^-factorization (or semihard) QCD approach, which is more preferable for 
small region than SPM. 

The /^-factorization QCD approach is based on Balitsky, Fadin, Kuraev, Lipa- 
tov (BFKL) [9] evolution equations. The resummation of the terms ln n (// 2 /AQ CD ), 

ln n (/i 2 /AQ CD ) ln n (l/:r) and a r g ln™(l/x) in the /c T -factorization approach leads to the 
unintegrated (depends from q T ) gluon distribution $(x, q^,/i 2 ) which determine the proba- 
bility to find a gluon carrying the longitudinal momentum fraction x and transverse momen- 
tum at probing scale fi 2 . The unintegrated gluon distribution reduce to the conventional 
gluon distribution xG(x,fj, 2 ) once the qr dependence is integrated out: 



To calculate the cross section of a physical process the unintegrated gluon distributions 
have to be convoluted with off mass shell matrix elements corresponding to the relevant par- 
tonic subprocesses [6-8]. In the off mass shell matrix element the virtual gluon polarization 
tensor is taken in the BFKL form [5]: 



Nowadays, the significance of the fey-factorization QCD approach becomes more and more 
commonly recognized [10]. It is already was used for the description of a wide class heavy 
quark and quarkonium production processes [6, 11-25], in particular, for the description 
heavy quark [6, 11, 14, 16, 20-23] and quarkonium [13, 15, 19, 24, 25] photo and leptopro- 
duction. It is notable that calculations in /^-factorization QCD approach results to some 
effects which are absent in other approaches. It is a more fast growth of total cross sections 
in comparison with SPM, a broadering of the pr spectra due to extra transverse momentum 
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of the colliding partons and other polarization properties of final particles in comparison 
with SPM. 

We point out that heavy quark and quarkonium cross section calculations within the SPM 
in the fixed order of pQCD have some problems. For example, the very large discrepancy (by 
more than an order of magnitude) [26, 27] between the pQCD predictions for hadroproduc- 
tion J /ty and T mesons and experimental data at Tevatron was found. This fact has resulted 
in intensive theoretical investigations of such processes. In particular, it was required to en- 
ter additional transition mechanism from cc-pair to the J /iff mesons, so-called the colour 
octet (CO) model [28], where cc-pair is produced in the color octet state and transforms into 
final colour singlet (CS) state by help soft gluon radiation. The CO model was supposed 
to be applicable to heavy quarkonium hadro and electroproduction processes. However, the 
contributions from the CO mechanism to the J /iff meson photoproduction contradict the 
experimental data at HERA for ^-distribution [29-32]. 

Another difficulty of the CO model is J /iff polarization properties in pp-inter act ions at 
Tevatron. In the framework of the CO model J /iff mesons should be transverse polarized at 
the large transverse momenta p^. However, this fact contradicts the experimental data too. 

The CO model has been applied earlier in analysis of the J /iff inelastic production exper- 
imental data at HERA [33, 34]. However, this results are not agree with each other [34]. It 
is notable that results obtained within the usual collinear approach and CS model [35-38] 
underestimate experimental data by factor about 2. 

Recently the first attempts to investigate a J/f polarization problem in pp-interactions 
at Tevatron in /^-factorization approach were made [24, 39, 40]. Also the theoretical pre- 
diction within semihard QCD approach [19] are stimulated the experimental analysis J /iff 
polarization properties at HERA conditions. However the further theoretical and experi- 
mental studies of this problem are necessary. 

Based on the above mentioned results we consider here the deep inelastic J /iff production 
at HERA in the framework of the CS model and the /cy-factorization approach with emphasis 
of a role of a proton gluon distribution function. Using the CS model and formalism of the 
projection operator [41], we investigate the Q 2 -, pf-, z-, y*- and TV-dependences of J /if! 
production on different forms of the unintegrated gluon distribution. Special attention is 
given to the unintegrated gluon distributions obtained from BFKL evolution equation which 
has been applied earlier in our previous papers [15-18]. For studying J /if! meson polarization 
properties we calculate the p^- and Q 2 -dependences of the spin aligment parameter a. 

The outline of this paper is as follows. In Section 2 we give the formulas for the inelastic 
J /iff electroproduction differential cross section in the /c^-factorization QCD approach. In 
Section 3 we describe the unintegrated gluon distributions which we use in our calculations. 
In Section 4 we give the formulas for J /iff electroproduction off mass shell matrix elements 
in the CS model. In Section 5 we present the results of our calculations. Finally, in Section 
6, we give some conclusions. 



2 Deep inelastic J/\I/ electroproduction cross section in the Ut- 
factorization QCD approach 

In this section we calculate total and differential cross section for inelastic J /iff electropro- 
duction via the photon-gluon fusion QCD subprocess (Fig. 1) in the /c T -factorization QCD 
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approach. Let us define Sudakov variables of the process ep — > e' J /if? X: 

Pj/y = aip e + f3ip p + PJ/M-T, P 3 = a 2 p e + foPp + PgT, , . 

qi = xip e + q 1T , q2 = x 2 p p + q 2 T, 

where 

2 2 2 n 2 2 -o2 2 2 /jN 

Pj/* = m J/*i Pg = °, ?1 = QlT = -Q > ?2 = ?2T> ( 4 ) 

Pj/* and p s are 4-momenta of the J/\I/ meson and final real gluon, qi and g 2 ar e 4-momenta 
of the initial virtual photon and gluon, Pj/^t, PgT, qir, qiT are transverse 4-momenta of 
these. In the center-of-mass frame of colliding particles we can write: 

p e = (E, 0, 0, E), p p = (E, 0, 0, —E), (5) 

where 

E =^ PI=P 2 P = Q, (Pe'flp) = |. (6) 

Sudakov variables are expressed as follows: 

m J/*T / s PgT i \ 

V s V s (7) 

a m J/*T ( v fl P S T ( v V ' 

Pi = — W- exp(-j/j/*), = —7=- exp(-y g ), 

where mj^ T = m 2 /^ + p 2 /^ T , mj/* is J/^ meson mass, yj/q and y 9 are cm. rapidities 
of J /ty meson and final gluon respectively. The inelastic J /if? electroproduction differential 
cross section has the following form: 

d<r(ep^e'J/VX) = — $(x 2 , q 2 T , /i 2 ) ^ dq 2 T da(e g* -> e' J/* </), (8) 

X 2 27T 

where 2 is initial BFKL gluon azimuthal angle, $(x 2 , q| T , Z^ 2 ) is an unintegrated gluon 
distribution in the proton, and 

dfr( e 0. _> e'J/if?g') = ^ £ l^lWe/ - e> J/V </) x 

(9) 

dhi d^pj/y d s p g (4) , 
X (2^2 (£5^ (2^ * fe + * " " P '"* " 

where J] I^IIhaI* 3 ^* — * e' J/if? g') is the off mass shell matrix element. In (9) J2 indicates 
an averaging over initial particles polarization and a sum over polarization of final ones. 

From (8) and (9) we obtain the following formula for the inelastic J /if? electroproduction 
differential cross section in the /^-factorization QCD approach: 



d.(ep ✓ j/» x) = — ) *'/• x 

x El^llHA(e«--e'J/* 9 0dp5 / „ I .dQ J d4r^^^, 

where 



(10) 



2tt 2tt 2vr 

xi = ai + a 2 , x 2 = (3 1 +(3 2 , qiT + q 2 T = Pj/*t + P s t, (H) 
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0i and are initial virtual photon and outgoing J/\P meson azimuthal angles, z = 

(pj/v " Vp)/{.1i ' Vp)- If we take the limit q| T — > and if we average (10) over the transverse 
directions of the vector q 2T , we obtain the formula for the differential cross section in the 
SPM: 

da(e P - e' J/* X) = ^ ^ ^ L _ ^ -^L_ * 2 G(* 2 , „ 2 ) x 

where X) |M|p M (e g — > e' J/\I/ (?') is the matrix element in the SPM. Here indicates an 
averaging over initial particles polarization and a sum over polarization of final ones. We 
average over the transverse directions of the vector q 2 T using the following expression: 

/ /^ $(X2 ' ^ £I M Isha(^* -e' J/*<f) = 
= x 2 G(x 2 , /i 2 ) E \M\ 2 PM (eg - e' J/tt</), 

where 

Z7r q 2T 2 



3 Unintegrated gluon distributions 

Various parametrizations of the unintegrated gluon distribution used in our calculations 
are discussed below. 

We used the results of a BFKL-like parametrization of unintegrated gluon distribution (so 
called JB parametrization) according to the prescription given in [44]. The proposed method 
lies upon a straightforward perturbative solution of the BFKL equation where the collinear 
gluon density xG(x,fi 2 ) from the standard GRV set [45] is used as the boundary condition 
in the integral form (1). The unintegrated gluon density is calculated as a convolution of 
collinear gluon distribution xG(x,fi 2 ) with universal weight factors: 



d V , (15) 



where 



as 
as 



J (2y/a s ]n(l/ V ) hnV/ql)), if c& < /i 2 
2y/asMVv)M<&/n 2 )), ifq^>/i 2 



(16) 



where J and J stand for Bessel function of real and imaginary arguments respectively, and 
as = Sas/n- The parameter as is connected with the Pomeron trajectory intercept: A = 
4aslii.2 in the LO and A = 4cSs'ln2 — Na s in the NLO approximations, where N ~ 18 [46, 
47]. The latter value of A have dramatic consequences for high energy phenomenology, 
however, some resummation procedures proposed in the last years lead to positive value of 
A ~ 0.2 — 0.3 [47, 48]. Therefore in our calculations with (17) we used only the solution 
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of LO BFKL equation and considered A as a free parameter varying it from 0.166 to 0.53. 
Pomeron intercept parameter A = 0.35 was obtained from the description of pt spectrum 
of D* meson electroproduction at HERA [20]. We used this value of the parameter A in 
present paper. 

Then we use another set (the KMS parametrization) [49] is obtained from a unified BFKL 
and DGLAP description of F 2 data and includes the so called consistency constraint [50]. 
The consistency constraint introduce a large correction to the LO BFKL equation: about 
70 of the full NLO corrections to the BFKL exponent A are effectively included in this 
constraint, as is declared in [51]. 

Finally, the third unintegrated gluon function used here is one proposed by Golec-Biernat 
and Wusthoff (GBW) which takes into account gluon saturation effects and has been applied 
earlier in analysis of the inclusive and diffractive ep-scattering data [52, 53]. It has the 
following form: 

$(x,q|,/i 2 ) = p^ — Rl( x ) q l e W (-R 2 (x)cfr), (17) 



Air 2 as 



where 



1 / rr \ A/2 



W = G^UJ ■ (18) 

and a = 23.03 mb, A = 0.288, x = 3.04 • 10~ 4 . 

In expression (1) the accounting of the contribution of the unintegrated gluon distribution 
to integral at low qf, region < qf~ < Qq was done by the collinear gluon density xG(x, Qq), 
where Qq = 1 GeV 2 for all gluon distributions. 



4 Off mass shell matrix element for eg* — > e' J/ty g' process 

In this section we give the formulas for J/\I/ electroproduction off mass shell matrix 
elements in the CS model. There are six Feynman diagrams (Fig. 2) which describe partonic 
process 7 g* — > J/\& g' at the leading order in as and a. In the framework of the CS 
model and nonrelativistic approximation the production of J/^f meson is considered as the 
production of a quark- ant iquark system in the colour singlet state with orbital momentum 
L = and spin momentum S = 1. The binding energy and relative momentum of quarks 
in the J/\I/ are neglected. Such a way mj/q = 2m c , where m c is charm mass. Taking into 
account the formalism of the projection operator [41] the amplitude of the process 7$* — > 
J/ty g' may be obtained from the amplitude of the process 7 (?* — > ccg' after replacement: 

^(0) 1 
v{Pc) u(p c ) -> J(pj/*) = — j= e(pj/y) (pj/vj, + raj/*) (19) 

2yjmj/^ ^3 

where i(pj/y) = e^(pj/^) 7^, e(pj/*) is a 4-vector of the J/^f polarization, 1/V3 is the color 
factor, ^(0) is the nonrelativistic meson wave function at the origin. In according with the 
standard Feynman rules for QCD diagrams we have: 



xSp 



M = e c g 2 e M (gi) e a (q 2 ) e p (p g ) x 
i, s u Pc-qi + m c - /v • p., ■ in,. 1 ( 20 ) 



(p c - qi) 2 - m 2 c (-Pc-PgY 
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+ 5 permutations of all gauge bosons. Here e M (gi), e^fe) are polarization vectors of the 
initial photon and gluon respectively, e M (p s ) is a 4-vector of the final gluon polarization. 
The summation on the J/ty meson and final gluon polarizations is carried out by covariant 
formulas: 

T,*<pj/*W(pj/*) = -r + (21) 

*Ze»(p g y(p g ) = -g<", (22) 

The initial BFKL gluon polarization tensor is taken in form (2). We substitute here the 
full lepton tensor for the photon polarization tensor (including also the photon propagator 
factor and photon-lepton coupling): 

E^teK^i) = 2^ + ■ (23) 

For studing J/ty polarized production we introduce the 4-vector of the longitudinal polar- 
ization e1(pj/q,) as follows [54]: 

= , {Pm ' Pp) { - t mj >* V \\ . (24) 

yj(pj/* -Pp) 2 - s \ m J/* (PJ/v ■ Pp) J 

In the limit s 3> fr?jm the polarization 4-vector satisfies usual conditions: (ei = —1 and 
(cl-Pj/*) = 0. The calculation of J2 I^Isha^* ->e! J/Vcf) and £ \M\l M (eg -> e' J/V g') 
was done analitically by REDUCE system. 



5 Numerical results and discussion 

In this section we present the theoretical results in comparison with available experimental 
data taken by the HI [42] and ZEUS [43] collaborations at HERA. There are three parameters 
which determine the common normalization factor of the cross section under consideration: 
J meson wave function at the origin ^(0), charmed quark mass m c and factorization scale 
jj,. The value of the J/^> meson wave function at the origin may be calculated in a potential 
model or obtained from the well known experimental decay width T(J/^/ — ► /i + /i~). In 
our calculation we used the following choice |\1/(0)| 2 = 0.0876 GeV 3 as in [55]. Concerning 
a charmed quark mass, the situation is not clear: on the one hand, in the nonrelativistic 
approximation one has m c = irtj/^/2 = 1.55 GeV, but on the other hand there are many 
examples when smaller value of a charm mass is used, for example, m c = 1.4 GeV [34, 56, 
57]. In the present paper we used both of charm mass values. Also the most significant 
theoretical uncertanties refer to the choice of the factorization scale /zp and renormalization 
one hr. One of them is related to the evolution of the gluon distributions $(x, qr,/ip), the 
other is responsible for strong coupling constant as (/^r)- As is often done in literature for 
simplicity we set them equal: fip = \ir — \i. In the present paper we used the following 
choice /1 2 = q| T as in [58] .0 

The calculations of the inelastic J/ty production cross section in the fc^-factorization 
approach have been made according to (10) and in the SPM — according to (12). The 

3 )The our preliminary results at [59] were obtained with p? = m 2 ,,^ + P 2 j/^, T - 
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p 2 j/q,T integration limits are taken as given by inelastic J/^f process requirements: 1 GeV 2 < 
Pj/*t < s/4 — mj^. According to (1) the q^ T integration region was divided into two parts: 
< qIt — Qo an d qIt — Qo- The calculations are made according to (10) for q^ T > Ql and 
for < c£ T < Ql we use £ |M|| M (e 5 -> e' J/* instead of £ I^IshaO e ' 9') 
(see also (12)). The choice of the critical value of parameter Q% = 1 GeV 2 is determined by 
the requirement that the value of as (q|r) be small in the region q^ r > lGeV 2 , where in 
fact a 5 (q 2 T ) < 0.26. 

The Q 2 , yj/q and z integration limits are taken as given by kinematical conditions of HI 
experimental data [42]. The first kinematical region is 2 GeV 2 < Q 2 < 80 GeV 2 , 40 GeV < 
W < 180 GeV, z > 0.2, M x > 10 GeV and the second one is p 2 7/ ^ r > 4 GeV 2 , 4 GeV 2 < 
Q 2 < 80 GeV 2 , 40 GeV < W < 180 GeV, z > 0.2, M x > 10 GeV. 

The results of our calculations are shown in Fig. 3 — 8. Fig. 3 show the Q 2 -, p 2 /^, T -, 
z-, yj/%- and iy-distributions of the inelastic J/\& meson production obtained in the first 
kinematical region at yfs = 314 GeV, m c = 1.55 GeV and Aqcd — 250 MeV. Curve 1 
corresponds to the SPM calculations at the leading order approximation with GRV gluon 
density, curves 2, 3 and 4 correspond to the /^-factorization results with JB, KMS and GBW 
unintegrated gluon distributions at Qq = 1 GeV 2 . One can see that results obtained in the 
CS model with fc-r-factorization (curves 2 — 4) agree in shape but underestimate the HI 
experimental data. The SPM calculation (curve 1) are lower the data by a factor 2. 

Fig. 4 show the z-, y}/^- and ^-distributions of the inelastic J/^ meson production 
obtained in the second kinematical region at \fs = 314 GeV, m c = 1.55 GeV and Aqcd — 
250 MeV. Curves 1 — 4 are the same as in Fig. 3. One can see that in this kinematical 
region the /cr-factorization approach (curves 2 — 4) describe the data better than in the 
first one (without cut p 2 /^ T > 4 GeV 2 ). 

Fig. 5, 6 shows the Q 2 -, p 2 ^ T -, z-, y*jiy- and VT-distributions of the inelastic J/ty meson 
production at ^fs = 314 GeV, m c = 1.4 GeV and Aqcd = 250 MeV obtained in the first 
kinematical region (Fig. 5) and in the second one (Fig. 6). Curves 1 — 4 are the same as 
in Fig. 3. One can see that a shift the charm quark mass down to m c = 1.4 GeV increases 
the cross section by a factor 1.5, and results obtained in the CS model with /c^-factorization 
(curves 2 — 4) agree with HI experimental data. 

The ^-distributions are described at z > 0.4 only. It is because other J/ty production 
mechanisms (such as resolved photon and/or colour octet contributions) may be impotant 
at z < 0.4 [60]. 

It is notable that the saturation unintegrated gluon distribution (GBW parametrization) 
does not contradict the HERA experimental data. 

Fig. 3 — 6 show that /cy-factorization results with m c = 1.4 GeV (in contrast with SPM 
ones) for inelastic J/ty electroproduction at HERA reasonably agree with HI experimental 
data. Hovewer, at m c = 1.55 GeV some additional transition mechanism from cc-pair to the 
J/ty mesons (such as CO model) need for full description of these data. 

We analyzed specially the influence of charm quark mass, m c , on the our theoretical 
results. We found that the main effect of change of the charm quark mass connects with 
final phase space of J/\l/ meson. In the subprocess matrix elements this effect is neglectable. 
However the value of m c = 1.4 GeV corresponds to the unphysical phase space of J/m state.0 

For studying J/ty meson polarization properties we calculate the p 2 -- and Q 2 -dependences 

4 )\Ve thank S. Baranov for the suggestion to study this problem. 
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of the spin aligment parameter a [19, 61] in the /c^-factorization approach and SPM: 



. . dalduj — 3 darjduj 

a(u) = -7^73 , / , , 25 

v ' da/du + da L /du v ' 

where <7l is the production cross section for the longitudinally polarized J/ty mesons, uu = 
Pj/*T' Q 2 ■ The parameter a controls the angular distribution for leptons in the decay 
J/\& — > fi + fi~ (in the J/^> meson rest frame): 



dcos8 



~l + acos 2 #. (26) 



Fig. 7 shows the parameter ol(pj/\&t)i which is calculated in the region 0.4 < z < 0.9 (a) 
and in the region 0.4 < z < 1 (b) at Vs = 314 GeV, m c = 1.4 GeV and A QCD = 250 MeV 
in comparison with preliminary experimental data taken by the ZEUS [43] collaboration at 
HERA. Curve 1 corresponds to the SPM calculations at the leading order approximation 
with GRV gluon density, curve 3 corresponds to the /c^-factorization QCD calculations with 
KMS unintegrated gluon distribution at Q\ — 1 GeV 2 . 

We note that it is impossible to make of exact conclusions about a BFKL gluon contri- 
bution to the polarized J/\& production cross section because of large incertainties in the 
experimental data and large additional contribution from initial longitudinal polarization 
of virtual photons. However at low Q 2 < 1 GeV 2 such contributions are neglected. This 
fact should result in observable spin effects of final J/^f mesons. As example, we have per- 
formed calculations for the spin parameter a as a function p 2 /^ T at fixed values of Q 2 for 
40 GeV < W < 180 GeV, z > 0.2, M x > 10 GeV at = 314 GeV, m c = 1.4 GeV and 
A QCD = 250MeV. 

The results of our calculations at fixed values of Q 2 = 0.1 GeV 2 , Q 2 = 1 GeV 2 , Q 2 = 
5 GeV 2 and Q 2 = 10 GeV 2 are shown in Fig. 8. Curves 1 and 3 are the same as in Fig. 7. We 
have large difference between predictions of the leading order of SPM and the /c T -factorization 
approach at low Q 2 < 1 GeV 2 , as it is seen in Fig. 8. Therefore if account of the NLO 
corrections (which is still unknown) will not change predictions of the LO SPM for a(p 2 ^ T ) 

at low Q 2 < 1 GeV 2 , experimental measurement of polarization properties of the J/\& mesons 
will be an additional test of BFKL gluon dynamics. 



6 Conclusions 

In this paper we considered J/^f meson production in ep deep inelastic scattering in the 
colour singlet model using the standart parton model at the leading order in as and the 
&T-factorization QCD approach. We investigated the z-, Q 2 -, p 2 --, y*-, H^-dependences of 
inelastic J/\I/ production on different forms of the unintegrated gluon distribution. The p 2 -- 
and Q 2 -dependences of the spin aligment parameter a presented also. We compared the the- 
oretical results with available experimental data taken by the HI and ZEUS collaborations 
at HERA. We have found that the /c T -factorization results with JB, KMS and GBW uninte- 
grated gluon distributions reasonably agree with HI experimental data only at m c = 1.4 GeV, 
|^(0)| 2 = 0.0876 GeV 3 and A QC d = 250 MeV. At m c = 1.55 GeV some additional transition 
mechanism from cc-pair to the J/ty mesons (such as CO model) need for full description 
of these data. The saturation unintegrated gluon distribution (GBW parametrization) does 
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not contradict the existing experimental data. The experimental study of a polarization of 
J/ty meson at low Q 2 < 1 GeV 2 should be additianal test of BFKL gluon dynamics. 

The authors would like to thank S. Baranov for encouraging interest and useful discus- 
sions. A.L. thanks also V. Saleev for the help on the initial stage of work. The study was 
supported in part by RFBR grant 02-02-17513. 
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Fig. 1. Diagram for ep — > e' J X process. 

Fig. 2. Feynmans diagram used for description partonic process 75 — > J/^f g' process. 

Fig. 3. The Q 2 -, p 2 /^ T -, -2-, and W^-distributions of the inelastic J/\I/ production 

obtained in the kinematical region 2 GeV 2 < Q 2 < 80GeV 2 , 40 GeV < W < 180 GeV, 
z > 0.2, M x > 10 GeV at ^ = 314 GeV, m c = 1.55 GeV and A QCD = 250 MeV. Curve 1 
corresponds to the SPM calculations at the leading order approximation with GRV gluon 
density, curves 2, 3 and 4 correspond to the /^-factorization QCD calculations with JB, 
KMS and GBW unintegrated gluon distribution at Ql = 1 GeV 2 . 

Fig. 4. The z-, y*ji^- and ^-distributions of the inelastic J/\I/ production obtained in the 
kinematical region p 2 J/<!)T > 4 GeV 2 , 4 GeV 2 < Q 2 < 80 GeV 2 , 40 GeV < W < 180 GeV, 
z > 0.2, M x > 10 GeV at yfi = 314 GeV, m c = 1.55 GeV and A QCD = 250 MeV. Curves 1 
— 4 are the same as in Fig. 3. 

Fig. 5. The Q 2 -, p 2 /^-, z-, y}/^- and W^-distributions of the inelastic J/ty production 
obtained in the kinematical region 2 GeV 2 < Q 2 < 80 GeV 2 , 40 GeV < W < 180 GeV, 
z > 0.2, M x > 10 GeV at = 314 GeV, m c = 1.4 GeV and A QCD = 250 MeV. Curves 1 - 
4 are the same as in Fig. 3. 

Fig. 6. The z-, y}/^- and W^-distributions of the inelastic J/ty production obtained in the 
kinematical region p 2 ^ r > 4 GeV 2 , 4 GeV 2 < Q 2 < 80 GeV 2 , 40 GeV < W < 180 GeV, 
z > 0.2, M x > 10 GeV at ^/s = 314 GeV, m c = 1.4 GeV and A QCD = 250 MeV. Curves 1 - 
4 are the same as in Fig. 3. 

Fig. 7. The parameter q;(pj/*t), which is calculated in the region 0.4 < z < 0.9 (a) and 
in the region 0.4 < z < 1 (b) at = 314 GeV, m c = 1.4 GeV and A QCD = 250 MeV. 
Curve 1 corresponds to the SPM calculations at the leading order approximation with GRV 
gluon density, curve 3 corresponds to the /c r -factorization QCD calculations with KMS 
unintegrated gluon distribution at Ql — 1 GeV 2 . 

Fig. 8. The spin aligment parameter a(p 2 ^ T ) at fixed values of Q 2 for 40 GeV < W < 
180 GeV, z > 0.2, M x > 10 GeV at = 314 GeV, m c = 1.4 GeV and A QCD = 250 MeV. 
Curves 1 and 3 are the same as in Fig. 7. 



